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ABSTRACT 


Investigation of oceanographic data acquired in the Columbia River 
estuary, which lies between the states of Washington and Oregon, 
revealed significant variation in sound velocity with respect to time and 
location. National Ocean Service (NOS) standards for echo-sounder | 
acquired soundings require the knowledge of sound velocity to within 
plus or minus 4 meters per second. To meet this requirement in the 
Columbia River estuary, methods were devised to allow prediction of 
sound velocity based on location and the height of tide. 

Confidence intervals associated with sound velocity predictions 
provided by regression analyses showed a substantial improvement in 
accuracy when compared to using a single average velocity for each 
location. Nevertheless, there were still some cases where the NOS 


requirement could not be met. 
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l. INTRODUCTION 


A. NOS HYDROGRAPHY 


“Hydrography is that branch of physical oceanography sles RK with 
the measurement and definition of the configuration of the bottoms and 
adjacent land areas of oceans, lakes, Biers. harbors, and other water 
forms on earth" [Ref. 1: p. 1-3]. The National Ocean Service (NOS) 
is the federal agency responsible for producing accurate nautical charts 
of nearshore areas of the United States and its territorial possesions. 
The principal objective of most hydrographic surveys conducted by NOS 
is to obtain basic data for the compilation of nautical charts with 
emphasis on the features that may affect safe navigation [Ref. 1: p. 
1-3]. For this thesis data analysis is limited to the study of sound 
velocity variations and the correction of echo soundings for NOS 
hydrographic surveys. The type of data used for this investigation is 
often the basis for other physical studies such as the dynamics of 
oceanic mixing and salt transport. No attempt will be made to analyze 


or explain such processes. 


1. Sound Velocity and Sound Speed 


Manuals and papers about hydrographic surveying often use 
the term "sound velocity" to designate what a physicist would define as 
“sound speed." The term “sound velocity’ will be used in this thesis 


as many of the references also use the term. 
2. Echo Sounding 


Echo sounders are the basic tool of the NOS hydrographer for 
determination of depth. Echo sounders measure the time required for a 
sound wave to travel from its point of origin to the bottom and return. 
They also convert this measured time to distance or depth [Ref. 1: p. 


AG-1]. The basic formula for echo sounder depth determination is: 


Zo = Vane 2 Cia) 


where Z is the depth, V is the calibrated velocity of sound for the echo 
sounder, and T is the time between transmission and reception of the 
sound signal. For depths to be accurate precise measurements of time 
and a knowledge of the mean, or harmonic, sound velocity (Chapter !1) 


are mandatory. 


B. THE VELOCITY OF SOUND IN SEA WATER 


Wilson [Ref. 2: p. 641] developed empirical equations based on 
laboratory measurements of sound velocity under tightly controlled 


conditions. The basic equation is: 
Vo 144 a en Pee cue cata (2) 


where V is the velocity of sound in meters per second, va is a 

~ correction for temperature, Vp is a correction fOr PRCSSUGe ys calome 
correction for salinity, and Vce7p is a simultaneous correction for the 
combined effects of temperature, salinity, and depth [Ref. 1: p. 4-75]. 
Wilson's equations (Chapter I!), which are currently used by NOS, 


were chosen to analyze the data for this thesis. 


C. NOS ACCURACY REQUIREMENTS 


The sound velocity at which NOS echo sounders are calibrated is 
800 fathoms per second. To obtain accurate depths for nautical charts 
the hydrographer must correct all soundings where the velocity of 
sound deviates from the calibrated value. There are several methods 
available for determining corrections for soundings, and they will be 
discussed in the following section. The NOS accuracy requirements 
[Ref. 1: p. 4-70] are: 


For use in correcting echo Sern ae the velocity of sound must 
be known with sufficient accuracy to ensure that no sounding will 
be in error by_as much as 0.25 percent of the depth from this 
cause alone. Therefore, the mean velocity must be Known to 
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within plus or. minus 4 meters per second. To satisfy these 
accuracy requirements, one must know the mean temperature of the 
water to an accuracy of plus or minus 1 degree Celsius and the 
salinity to within plus or minus 1 part per thousand. A sufficient 
number of salinity and temperature observations must be made so 
that the velocity can be determined within the specified accuracy 
over the entire area sounded. 


D. PROCEDURES FOR CORRECTING SOUNDINGS 


Echo soundings can be corrected by either a direct comparison or 
an oceanographic procedure. Direct comparisons involve comparing echo 
soundings to known depths; this does not require a Knowledge of the 
sound velocity. Procedures which involve the calculation of sound 
velocity from measurements of water characteristics are classified as 
“oceanographic and limnologic determinations” [Ref. 1: p. 4-74]. For 
this thesis the direct measurement of sound velocity by a sound 


velocimeter is also classified as an oceanographic procedure. 
1. Direct Comparison Systems 


Two procedures used to determine depth corrections by the 
direct comparison method are the bar check and the vertical cast. In 
these the effects of sound velocity variation and instrument error are 
included in the results [Ref. 1: p. 4-71]. 


a. Bar Check 


The bar check procedure involves lowering a bar over 
the side of a survey vessel and positioning the bar under the 
transducer at a Known depth indicated by markings on the lines 
supporting the bar. Comparisons between the depth as indicated by 
the echo sounder and the Known depth are made at 5- or 6-foot 
increments as the bar is lowered and raised though the water column. 
Due to the handling of the bar, this procedure is usually limited to 
smaller vessels and the quality of the results is a function of sea 
-conditions. Observations during rough seas are subject to unacceptable 


magnitudes of error [Ref. 1: p. 4-71]. 


b. Vertical Cast 


The vertical cast requires simultaneous sounding of the 
bottom using an echo sounder and a lead line. This procedure, which 
can be accomplished by practically all survey vessels, is not as 
accurate as the bar check, but vertical casts can be used to determine 
or verify instrument errors. Where sounding sheets junction, there are 
often discrepancies and soundings fail to agree. Vertical casts can 


sometimes help resolve such problems [Ref. 1: p. 4-72]. 
2. Oceanographic Sound Velocity Determinations 


The more traditional oceanographic procedures involve sound 
velocity determination by measurement of the variables which enter 
Equation 1.2, namely temperature and salinity as functions of depth. 
The Nansen cast and electronic sensors such as the CTD (conductivity, 
temperature, and depth) or the STD (salinity, temperature, and depth) 
are examples of equipment used for these procedures. An alternative is 
the direct measurement of sdund velocity using a sound velocimeter. 
The electronic sensors record data almost continuously with depth. 

[Ref. 1: p. 4-74]. 

To obtain velocity corrections with Nansen, STD, or CTD 
data, Umbach [Ref. 1: p. 4-75] recommends using the "summation of 
layers’ method, which extends from the concept of mean or harmonic 
sound velocity. 

Most of the electronic sensors which are currently in use yield 
precise measurements in excess of NOS accuracy requirements for sound 
velocity corrections if properly calibrated and maintained [Ref. 1: p. 
4-74]. According to NOS standards all electronic sensors must be 
periodically checked for accuracy against Nansen cast or equivalent 


observations. 


E. SOUND VELOCITY VARIATION 


The effect of changes in temperature, salinity, and pressure on 


sound velocity depends upon the initial conditions of these variables. 


Wilson's equations are used to illustrate in Table | the effect of 


changing variables by one unit. 


TABLE | 
Variation of Sound Velocity 


PRES TEMP SAL V delta V 
db ae g/kg m/s m/s 
Oe 15.00 30.00 1502.90 
eye Lal Gye) ee) Sele. 1499.30 3.60 
iL Ce 5.00 30.00 1464.27 
iL 010 HOO sO. 00 1460.03 4.25 
' 
1.00 10.00 35.00 1490.76 | 
iL, O16 ieee "34, 00 1489.48 1.28 
oo 1,0) (016, iL. 06 1449.01 
ele LC), (XO 0.00 1447.84 ieee 
60.00 LOu@Ou 30,00 1485.33 
59.00 iO OC. SO, 00 1485.32 0.02 
i @O 10400 30506 1484.38 
6700. 10.00 2000 1484.36 0.62 


Velocity (V) was computed using Wilson's equations. 


Corrections for .sound velocity variations in estuaries can present 
the hydrographer with a difficult task. According to Cameron and 
Pritchard [Ref. 3: p. 306], an estuary is a semi-enclosed coastal body 
of water having a free connection with the open sea and within which 
the sea water is measurably diluted with fresh water deriving from land 
drainage. Estuaries are extremely dynamic systems. Topography, 
river flow, and tidal action are important factors that influence the rate 
and extent of mixing salt and fresh water [Ref. 4: p. 4]. For the 


hydrographer, estuaries are areas where sound velocity varies not only 
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with location but also in time due to the tidal cycle and variations in 
river flow. Although no specific procedures are given, Umbach 


[Ref. 1: p. 4-75] notes the possible complexity of such areas: 


Under special circumstances such as . .. in areas of extensive. 
estuarine discharge of fresh water . . . series of regional velocity 
curves are required. Regional curves must be carefully studied to 
determine how they can be best grouped by area or time or to 3 
permit drawing of average curves. 


F.. OBJEGIEVES 


One objective of this thesis was to investigate the magnitude of 
sound velocity variation associated with significant changes in water 
characteristics within an estuary. 

A second objective was to investigate the feasibility of determining 
velocity corrections according to NOS standards in areas where extreme 
changes in water characteristics occur with respect. to space and time. 

A third objective was to write a FORTRAN program to enable 
further testing by the field hydrographer equipped with a personal 


computer. The use of commercial subroutines has been avoided. 


Il. METHODOLOGY AND PROCEDURE 


Sa AReEASOF STUDY 


The Columbia River estuary which discharges into the Pacific 
Ocean at latitude 46°1S'N, longitude 124°00'W has variations in water 
characteristics which make it an excellent area for investigation. 
According to Dyer [Ref. 4: p. 96], the Columbia River estuary is 
characterized by an unusual combination of large river discharge and 
strong tidal currents. Recent investigations have revealed differences 
in salinity of up to 30 grams per kilogram at individual stations during 
a single tidal cycle [Ref. 5: p. 104]. Due to this large difference and 
the availability of data from a recently completed survey, this region 


was chosen for investigation. 


B. DATA SOURCES 


To analyze the problem completely, water characteristics 
(temperature and salinity as functions of depth) must be known at 
stations spaced throughout the region, and each station occupied during 
a complete lunar tidal cycle. Variations caused by seasonal changes in 


fresh-water runoff should be examined thoroughly. 
tec Columbia River Estuary Data Development Program 


The Columbia River Estuary Data Development Program 
(CREDDP) was a 5-year, federally funded study which included 
approximately 2 years of data acquisition for investigation of physical 
and biological characteristics. CREDDP was designed to meet the needs 
of two major groups. One group consisted of local, state, and federal 
government agencies involved in planning and permitting activities. 

The second group included research scientists and educational 


institutions [Ref. 6: p. v]. 


iS 


One accomplishment of CREDDP was the acquisition of current 
velocity, conductivity, temperature, and depth (VCTD) data for a 
better understanding of physical processes such as currents, sediment 
transport, and salinity intrusion. VCTD profiling was carried out for 
13 consecutive days from 17 stations for periods of from 12 to 48 hours 
per station. Between October 16 and 27, 1980, approximately 1600 
VCTD casts were made. Of these, 458 casts were made available to 
users [Ref. 6: p. L-8]. 

Data were supplied to NPS by the U.S. Army Corps of 
Engineers (COE) at Portland, Oregon, on two 9-track magnetic tapes 
[Ref. 6: pp. 1-5]. 

For an analysis of sound velocity variation, VCTD casts from 
11 of the 17 stations were investigated. The stations, along with their 
positions and times are listed in Table Il. The plot of the geographical 
positions of the stations is shown in Figure 2.1. These stations 
covered the area where the maximum sound velocity variation would be 


expected as deduced from salinity intrusion diagrams of the estuary 
(Ref. 5: p. 104]. 


Dh NOS Bar Check Data 


NOS bar check data were used for a sound velocity variation 
analysis. The descriptive report for survey H-8422 listed the day each 
bar check was taken and the correction for each 5-foot increment of 
depth. Also given was a general term for the stage of the tide such as 
ebb, flood, high water, or low water. Data for 46 bar checks were 
acquired between May 9 and July 16, 1958 [Ref. 7: p. 1]. 


Sn NOS Predicted Tide Data 


NOS predicted tide data were obtained for correlation and 
regression with the CREDDP and the H-8422 bar check data. Predicted 
tides were used because field hydrographers rarely have the time or 


the proper tools for a detailed analysis of real tide data. 


TABLE II 
Columbia River Estuary Data Development Program 


STATION (STN) INFORMATION 


PCOUoln TON hi Mis er 
STN FOS TIION START EN 


2N _46°15'45"N 10/18/80 0930 10/19/80 0230 
Ware @ee sly EE ae! 

2S _46°15'13"N 10/18/80 1230 10/19/80 0930 
124°03'20"w Pais Bee) 

4NA 46°14'O03"N 10/20/80 2115 10/21/80 1745 
123°56'09"W 

4NB 46°14'47"N 10/20/80 2105 10/21/80 1730 
123°55'53"W 

4SA _46°12'42"N 10726/80 0007 10/27/80 0832 
123°57'09"w ee Chat 

4SB 46°12'S1"N 10/25/80 2350 10/27/80 0820 
123°56'S5"W LOL HONS 

5NA 46°13'55"N 10/20/80 1733 10/22/80 0930 
oe en 7 ay , | 

SNB 46°14'18"N 10/20/80 2035 10/22/80 0913 
e253 14 "w 

SNC 46°14'07"N 10/21/80 0935 10/22/80 1003 
123°53'14"W jas beset 

SSA 46°11'55"N 10/24/80 0830 10/25/80 2331 
123°48'38"Ww 

SSB 46°12'05"N 10/24/80 1126 10/25/80 1130 
123°48'54"W al Goa 


17 


COLUMBIA RIVER ESTUARY 
CREDDP STATIONS 


| WASHINGTON 
oe =) 
ees 





Figure 2.1 CREDODP Stations. 


C. CREDDP DATA ANALYSIS 
ie COE V CT OD Pataslanes 


COE VCTD data tapes contained pertinent information such as 
header data, time, temperature, depth, conductivity, and salinity. 
Mercier [Ref. 6: p. L-8] shows a breakdown of the data components for 
one record of data. The tapes contained more data records than were 
necessary for this investigation. The VCTD sampling device recorded 


data every 0.04 seconds, thereby making the average number of 
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records for a single cast approximately 1500. Having 270 casts to 
analyze, a method was devised to reduce the data to a more efficient 
size. A computer program was formulated to input the data for each 
cast and then output the first record and the records which 
corresponded to the nearest whole meter of depth. For example, a 
10-meter cast would theoretically have 11 records of data; the first 
record would represent the surface layer, and the remaining 10 records 
would be for 1-meter increments. In accordance with the NOS policy of 
using bar check data measured at 5- or 6-foot increments, data were 
eliminated except at Il-meter increments. A sample of a profile after 
editing is shown in Table III. Standard depths which represent depths 
of 1, 3, 5, 7, 9, 11, and 13 meters, were adopted to further simplify 
data processing and description. 

Several drawbacks were encountered using the CREDDP data. 
Because of the shallow depths of CREDDP stations, many of the casts 
were inadequate for correcting soundings in the deeper areas of the 
estuary. Some casts did not include surface layer data (0 to 2 


meters). 


ue Sound Velocity Computations 


a. Wilson's Equations 


Sound velocity was computed for each record by Wilson's 
equations which use temperature, pressure, and salinity for the 
computation. In the following equations V is sound velocity in meters 
per second; T is temperature in degrees Celsius, where -4 < T < 30 
°C; P is pressure in kilograms per square centimeter, where 1 < P < 
1000 kg/cm?’ and S is salinity in grams per kilogram, where 0 < §S < 37 
e/ kon 


VE Vn Ne * Ver (2.1) 
Deo 04 4532x10 “T+ - oy 


2.6045x10°3T3 + 7.9851x107°T4 


i 


1.60272x107 Pp + 1.0268x107>p2 + (2.3) 


V = 
Pp 
3.5216x10-9p3 - 3.3603x10"!2p4 
Vg = 1.39799(S-35) * 1.69202x10°3(S-35)? (2.4) 
Vetp = (S-35)(-1.1244x10°2T + 7.7711x10"/T + (295) 


7.7016x1079P - 1.2943x107/p2 + 
3.1580x10°8pT + 1.5790x10-9PT2) + 
P(-1.8607x107°4T + 7.4812x107°T2 + 
4.5283x10°8T3) + p2(-2.5294x107/T + 
1.8563x10°9T2) + p3(-1.9646x107 !9T) 
The standard deviation for the velocity is 0.30 meters per 
second [Ref. 8: p. 1357]. 


TABLE 111 
4SA--Cast 6--Sound Velocity Variables 
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b. Mean Sound Velocity 


Because echo sounders are time measuring devices, 
distance or depth is one-half the round trip travel time of a sound 
pulse multiplied by a suitable mean velocity. The mean sound velocity 
(MSV) is given by 


MSV = Z/T C235) 
where Z is the distance traveled in time T. To correct for sound 
velocity variations, the water column ts divided into a series of layers, 


each with thickness AZ, and an associated velocity V;. Therefore, the 


time interval for the sound wave to pass through the ith layer is: 

AT, = AZ,/V; | (2.7) 
Summing all layers results in: 

n 

1=1 | 
OG 

T = faz/v | (2.9) 

0 jun 0 
Substituting into Equation 2.6 results in: 

MSV = Z/({SdZ/V) (2.10) 
Equation 2.10 was evaluated using the trapezoidal method [Ref. 9: p. 
86]. NOS uses mean sound velocities in the “summation of layers” 
method described by Umbach [Ref. 1: p. 4-74] and Wallace [Ref. 10: 
p. RK530.1]. 


Two sound velocity listings computed using the data in 


Table II! are shown in Table IV. The sound velocity column was 
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computed using Equation 2.1, and the mean velocity column was 


computed using Equation 2.10. 


TABLE IV 
4ASA--CAST 6--Sound Velocity Data 


DATE: 26 OCTOBER 1980 
TIME: 0520 (PDT) 


SOUND MEAN 
DEPTH VEROCI Eye Sheer 
m m/s m/s 
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To simplify the analysis of the variation of mean sound 
velocity as a function of time, a computer program was written to 


arrange the data in individual arrays for each depth at each station. 
3. Mean Sound Velocity and Tide Height 


Inspection of the mean sound velocity arrays suggested that 
velocities were maximum at high tide and minimum at low tide. To test 
this, a correlation and regression analysis between the CREDDP derived 


mean sound velocities and predicted tide heights was performed. 


age 


a. Predicted Tide Heights 


To perform the correlation and regression analysis, a 
method had to be developed to calculate the height of the tide at any 
desired time. One method would be to obtain the harmonic constants 
for prediction of the Columbia River estuary tides from the NOS Tides 
and Water Levels Branch. These constants would then be used in the 
appropriate equations [Ref. 11: p. 123]. From a computer programming 
standpoint this method would be the simplest. A second method would 
be to formulate a computer program based on the graphical method 
(Figure 2.2) described in the NOS Hydrographic Manual [Ref. 1: p. 
4-68]. This method was chosen because normally the field 
hydrographer is not supplied with the harmonic constants for project 
areas. However, the tide tables are an essential tool of the 
hydrographer. 

The present NOS HYDROPLOT computer program for tide 
prediction lacks documentation and is not designed to solve this 
particular problem [Ref. 10: p. 500.1]. Therefore, a program was 
written. 

Points A and E (Figure 2.2) represent the predicted low 
and high tide, respectively. The straight line between these points 
was divided into four equal parts by points B, C, and D. Point B’ was 
plotted directly below B at one-tenth the total range, and point D was 
plotted directly above D using the same distance [Ref. 1: p. 4-68], 
forming three regions for the solution of predicted heights between 
points A and E. For the region between points A and B’', the curve 
was assumed to be parabolic, and in order to calculate the tide height 


for a designated time in this region, the equation 


Pee bi. ic (2.11) 


was used, where H is the tide height and T is the designated time. To 


calculate the coefficients a, b, and c; the equations 


Ha = aT 53 bT a * 3G (2<12) 


Sa Se male as 
NOS ANALYTICAL METHOD 
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Figure 2.2 Predicted Tide Graph. 
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Hp: = aTg + bTg: + c (22) 
0 =2aT, +b (2.14) 


were used. For times that lie between points D' and E, the same 
method was used. Between points B and D' the change in height was 
considered to be linear with respect to time and can be computed by 


the equation 
where m is the slope defined by 
m= Cec hlie:)/ Gina T a') (25 16) 


Using the program based on the above equations, tide heights were 
calculated for the times of all VCTD casts used in the analysis. "Data 
Array’ means a set of sound velocity and predicted tide height data 
that shows the variation of both variables with time for a particular 
station and depth (Table V). 


aye Data Plots 


Computer generated plots of predicted tide height and 
observed mean sound velocity (Figure 2.3) were made for all stations 


using mean velocities based on standard depths. 
om Correlation Analysis 


A correlation analysis routine for mean sound velocity 
and tide height was run for all data arrays. In addition to computing 
the correlation coefficient, the routine also computed the averages and 
the standard deviations of the mean sound velocity and the predicted 
tide height for each data array. Averages were computed by the 


equation 


n 
xX = D>. x-/n (247) 
ie 


ZS 


TABLE V 
4SB--Data Array--Depth = 7 m 
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where x; is the jth observation and n is the number of observations 
[Ref. 12: p. 124]. Standard deviations [Ref. 12: p. 128] were 


computed using the formula 
2 
s = 2, OG-®) /(n-1) (2.1) 
l= 
Computation of the simple correlation coefficient [Ref. 12: p. 258] of 


the mean sound velocity and predicted tide height arrays was done with 


equation 


(2.13) 
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Figure 2.3 


A sample of the results of a correlation analysis is listed in Table VI. 


TABLE VI 
4SB--Statistical Data--Depth = 7 m 


AVERAGE TIDE HEIGHT 
STANDARD DEVIATION OF TIDE HEIGHT 


AVERAGE MEAN VELOCITY 
STANDARD DEVIATION OF MEAN VEECGIT 


CORRELAT ION Cena ciEne 





d. Regression Analysis 


To determine the feasibility of predicting mean sound 
velocity, a first-order (or linear) regression analysis was performed on 


all data arrays. The data were fitted to the equation 
y=a+ bx (222.03) 


where y is the estimated mean velocity, a is the coefficient, b is the 
the intercept, and x is the tide height. The coefficient and intercept 


estimates were determined to minimize the sum of the residuals 
n n 
Bes 2 (¥-yp (2.21) 


where Y, is the observed mean sound velocity and yi is the estimacea 
velocity [Ref. 13: p. 533]. Determination of the coefficients and the 


intercepts was accomplished by the normal equations 


n n 
+ b = fe : 
ant b 2x > 2 Yi (2.22) 
n 2 n 
ex, 8 DO xX = 2» xy (223) 


where n is the number of observations for a particular station and 


depth [Ref. 13: p. 533]. The sum of the squares of the residuals is 


The standard error of estimate is 


see =/(SS / v) (2523) 


where v is the degrees of freedom given by 
v=ezp-n (2.26) 


where n is the number of sample observations and p is the number of 
estimated variables [Ref: 14: p. 131]. For z first-order analysis, there 
are two estimated variables. Results of a first-order analysis are listed 
m Table Vil. 

Confidence intervals associated with predicted mean 
velocities obtained by equations derived from the linear regression 
analysis were calculated to determine the accuracy of predictions. 
Confidence interval computation requires the standard error of estimate, 
the value of the student-t (t) distribution which is based on the 
desired percent of confidence, and the degrees of freedom for the 
sample. For this thesis, 90-percent confidence intervals were computed 


using the equation 


Omar oeyeeego)(l *(1/n) * C22) 


n 
1G sae 3, (x)-%)7] see 
T= 


Zo 


TABLE VII 
4SB--Linear Regression--Depth = 7 m 


COEFFICIENT = ZS 
< LNZERCEE — 1468.0 
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A RESIDUAL IS THE DIFFERENCE BETWEEN THE OBSERVED 
MEAN VELOCITY AND THE PREDICTED MEAN VELOCITY. 


where Cl is the confidence interval for a single estimated velocity (y;), 
x, is the tide height, X is the mean tide height, and n is the number of 
sample observations [Ref. 15: p. 30]. 
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Several plots of the data arrays indicated a curvilinear 
relationship between mean sound velocity and predicted tide height. 
Therefore, second- and third-order polynomial regression analyses were 
performed on all data arrays to provide a better fit of the data, 
thereby increasing the accuracy of the mean sound velocity prediction. 
The normal equation method used in the first-order analysis was 
extended for calculation of a second-order polynomial. To check the 
results of the programs for the normal equation method, a subroutine 
using the orthogonal polynomial method was obtained from the NPS 
Computer System Library [Ref. 16: p. 8]. A third-order regresssion 
analysis was also performed by the orthogonal polynomial subroutine to 
investigate the possibility of increasing the accuracy of prediction and 
fit. Based on results of the investigation and check (Chapter III), the 
first-order and second-order normal equation methods were designated 
as the primary methods of analysis. The normal equation method 
appeared to be the simplest to program. The second-order analysis 


uses the equation 
2 : z Z: 7.98 
y = a ay x ay X (2726) 


where y represents the estimated mean sound velocity and x is the 
predicted tide height [Ref. 13: p. 535]. Coefficients (ap, a), and as ) 
were determined to minimize the sum of residuals by the following 


normal equations 


ee ht 2B a (2.29) 
n n 
ag 2% 4 2 x. ay 2 x PS x, Y (230) 


n n n n 
Dee Sie 4 _ 2 
ag 2 ne 21.2, ne a, 2% 2 x Y (223),) 
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where n is the number of casts for a particular station and depth 
[Ref. 13: p. 536]. Results of a second-order analysis are listed in 
Table VIII. 


ou 


TABLE VIII 
4SB--Polynomial Regression--Depth = 7 m 
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Confidence intervals for predicted velocities were 
computed using Equation 2.27. The degrees of freedom for a 
second-order analysis was less than a first-order analysis because of 


the addition of one estimated variable. 
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D. NOS BAR CHECK DATA 


The NOS bar check data consisted of 46 bar checks which were 
acquired to determine the corrections for soundings of survey H-8422. 
The format of the data was a listing of corrections for 5-foot increments 
of depth. The initial analysis included determination of the maximum, 


minimum, and mean correction for each depth (Table IX). 


TABLE IX 
Bar Check Correction Data For Survey H-8422 
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For each depth of a bar check there were two corrections listed; 
the first correction was recorded as the bar was being lowered, and the 
second during raising. These two corrections were averaged to obtain 
a single value which was used for a correlation analysis with predicted 


tide height. 
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A plot of the bar check and tide height data (Figure 2.4) was 
made and visually examined. The plot indicates a poor correlation 
between bar check corrections and predicted tide heights. For 
verification, a correlation analysis (Table X) was performed using 


methods previously discussed. 


TABLE X 
Bar Check Statistical Data for Survey H-8422 


MEAN TIDE HEIGHT 
STANDARD DEVIATION OF TIDE HEIGHT 


MEAN CORRECTION 
STANDARD DEVIATIONS=G©® CORRECIIGT 


CORRE LAT TOP COPE ere Nt 
DEPTH OF. OBSERVATIONS 





Because of the poor correlation, actual tidal data were obtained to 
investigate the possiblilty that the results were related to the predicted 
.tides. No improvement in the results was accomplished by using actual - 
tide data. No further manipulation of the NOS bar check data was 
carried out because of the number of unknowns associated with the 
acquisition procedures and because the amount of speculation necessary 


could not lead to meaningful analysis. 
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Ili. DISCUSSION OF RESULTS 


A. CREDDP DATA 


Data acquisition procedures of CREDDP made an investigation of 
mean sound velocity variation possible, but there are factors which 
make the study incomplete. Because of this, a complete model for mean 
sound velocity prediction and correction based on methods presented in 
this thesis was not fully developed. However, the primary components 
of such a model were devised and will be explained in this chapter. 
The factors which make the investigation incomplete and 
reaommendations concerning verification of ideas and methods presented 
in this thesis will be addressed in Chapter IV. 

Investigation of the CREDDP data revealed that sound velocity in 
the Columbia River estuary varies with position, depth, and time. To 
analyze the problem the study area was divided into two primary 
sections, the north channel and the south channel (Table XI). The 
‘river mile’ distance represents the shortest deep-water path from buoy 
"T" at the river entrance to the station. The south chanmeleewhicomia 
maintained for commercial marine traffic, has a river mile scale which 
can be found on NOS chart 18521, whereas north channel stations were 
scaled using chart 18521. 

Two categories of mean sound velocity variation are considered. 
The term “horizontal variation’ designates differences in mean sound 
velocity between stations; i.e., how the velocity changes with location. 
The horizontal variation analysis is needed to develop zones or regions 
Where velocity correctors may be grouped. The term ‘temporal 
variation’ designates the change in mean sound velocity with time at 
each station. With temporal variation, mean sound velocity correctors 


for each zone change with time. 
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TABLE XI 
Channel Data 


north channel south channel 


Station River Mile Station River Mile 


2N 
4NB 
5NB 





1. Horizontal Variction 


For the analysis of horizontal variatian the two requirements 
made were: (1) that for each depth there be but one horizontal 
variation analysis and; (2) the average of all mean velocities obtained at 
the depth being analyzed be designated as the average velocity. Use 
of the average velocity allows each station to have a single velocity for 
each depth. The average velocity presents a problem if more 
observations were made on one end of the tidal spectrum than the 
other. For example, if more casts were acquired during periods of 
high water than periods of low water, the average velocity would be 
higher than if the majority of casts were acquired during low water 
periods. One of the CREDDP objectives was to acquire casts on an 
hourly basis. This objective was met in most cases, but there were 
some instances where the time interval between successive casts was 
greater than 1 hour. For example, station 4SB was chosen over station 
4SA for designation of the south channel because there were no casts 
acquired at station 4SA during the higher high tide. Another potential 
oroblem would be introduced if there were significant variations in the 
river flow. According to Jay [Ref. 5: p. 114], river flow remained low 
for the acquisition period. An example of horizontal variation of the 
average velocity is shown in Figure 3.1. These plots indicate the rate 


of horizontal variation of the average velocity is approximately the same 


of 


for the north and south channels. Data for the 11- and 13-meter 
depths were limited and did not cover the range of tide well in many 
cases. 

Regional zoning of mean sound velocity. correctors based solely 
on horizontal variation can be accomplished with little difficulty. Using 
the depth with maximum horizontal variation, zones can be drawn based 
on the NOS criteria of plus or minus 4 meters per second. A major 
problem still exists in that the temporal variation of mean sound velocity 
has not been accounted for when only the horizontal variation of the 


average velocity is considered. 
2. Temporal Variation 


An analysis of temporal variation of the mean sound velocity 
demonstrated that proper correction ‘of soundings in the Columbia River 
is a multi-dimensional problem. Sound velocity variables change 
significantly with time at the same location, and the magnitude of 
Variation can depend on a number of cyclic and noncyclic phenomena. 
Sound velocity variables change with the tide, and the magnitude of the 
variation is affected by monthly tidal. events such as spring and neap 
tides. The temporal and horizontal variations of velocity for the 
stations of the north and south channels are shown in Figure 3.1. 

Data for stations 2N and 2S were acquired during a neap tide 
period when the predicted tide range was minimum, whereas data 
acquisition for stations 4SA, 4SB, 5SA,and 5SB was accomplished 
during a spring tide period. Data for stations 4NA, 4NB, 5NA, 5NB, 
and SNC were acquired during a period of average tide range. As 
illustrated by Figure 3.1, stations 4NB and 4SB have about the same 
average velocity at the standard depths, but with respect to temporal 
variation, the range of velocity for station 4SB was significantly larger 
than that found at station 4NB. According to Jay [Ref. 5: pp 
114,117], the neap tide period was a one of minimum vertical mixing and 
maximum stratification, whereas during the spring tide period, maximum 
mixing, minimum stratification, and maximum salinity intrusion occurred. 


Examination of the casts for stations 4NB and 4SB (Table XII) 
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confirmed the previous statement. The cast for station 4NB acquired at 
predicted low tide showed a significant difference in temperature when 
compared to the low tide cast of station 4SB. The warmer temperatures 
of the low tide cast caused the the minimum velocity to be higher at 
station 4NB than at station 4SB. Comparison of the high tide casts 
indicated the surface layer salinity was 6.82 parts per thousand less at 
station 4NB than at station 4SB. The relatively large salinity 
difference for the high tide casts caused the maximum mean velocity to 
be lower at station 4NB when compared to the maximum mean velocity at 
station 4SB. Therefore, the increased energy associated with a larger 
tide range can increase the mean sound velocity range for certain areas 
in two ways. First, increased mixing associated with the larger tide 
range can moderate temperature for low tide cast, thereby causing 
minimum mean velocities to be lower. Second, increased salinity 
intrusion associated with a larger tide range can cause maximum mean 
velocities to be higher when compared to maximum mean velocities 


observed during an average or neap tide range. 
a. Correlation of Tide Height and Sound Velocity 


The effort to correct mean sound velocity for temporal 
variation began with a correlation analysis of predicted tide height and 
mean sound velocity. The results (Table XIII) indicate a strong 
relationship between tide height and mean sound velocity for most of 
the CREDDP data. Correlation in the 1l-meter layer was usually not as 
strong as the deeper layers for most of the stations. Dashed line 
entries in Table XIII indicate insufficient data for a correlation 


analysis. 
b. Regression Analysis 


First- (linear), second-, and third-order polynomial 
regression analyses were performed for each standard depth of the 
CREDDP sampling stations. The objective of these analyses was to 
provide a numerical method of predicting mean sound velocity using the 


predicted tide height. These methods would be considered completely 
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TABLE XII 
The Effect of Tide Range on Mean Sound Velocity Range 
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successful if the confidence interval for predictions based on all 
possible tide heights was less than, or equal to, plus or minus 4 meters 


per second. The employment of three different orders of polynomial 
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regression analysis was for comparison of fit, as well as comparison of 
the predicted velocities and associated confidence intervals. 

A potential source of error may exist when using 
predictions calculated from tide height values which are outside of the 
range of values used for the initial regression analysis. For example, 
if all casts were acquired during a neap tide period, such as was the 
case for stations 2N and 2S, velocities predicted by the regression 
analysis for high or low tide heights of a spring tide period are subject 
to error. For stations 4SB and 4SA, predictions were calculated using 
a 9.0-foot tide height; station 4SB (Figure 3.2} was occupied during 
higher high tide, whereas station 4SA (Figure 3.3) was not occupied 
during this period. The second- and third-order analysis predicted 
velocities which probably were not physically possible at station 4SA. 
Therefore, a basic objective in using polynomial regression models for 
mean sound velocity predictions and corrections should be the 
determination of the maximum and minimum velocities expected, and the 
minimum range of tide where these velocities could occur. The most 
accurate method for identification of these values would be to measure 


the sound velocity variables at the same location covering more phases 
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of the lunar cycle. The danger of predicting outside of the observed 
range is illustrated by Figure 3.3, whereas Figure 3.2 shows that 
predicting within the observed range can be done accurately. 

Determination of the most efficient polynomial regression 
analysis to use for prediction of mean sound velocity was accomplished 
by examination of plots of the data arrays and comparison of the 
calculation of predicted velocities and confidence intervals provided by 
each analysis. 

(1) Data Array Plots. Examination of the data array 
plots indicated that the first-order (linear) analysis would be sufficient 
for most of the data arrays, although some of the arrays did suggest a 
‘curvilinear trend. Curvature could be detected in the plots of stations 
4NA, 4SB (Figure 3.2), 4SA (Figure 3.3), and 4NB (Figure 3.4), with 
maximum curvature in the plots of stations 4SA and 4SB. This may be 
an indication that the degree of curvature is proportional to the range 
of tide since data for stations 4SA and 4SB were acquired during the 
spring tide period, whereas data acquisition for stations 4NA and 4NB 
was carried out during a period of average tide range. This 
hypothesis cannot be verified without further investigation. 

(2) Predicted Velocities. A computer program was 
written to predict the mean velocity for five different heights of tide 
with the maximum and minimum heights being the same as those used in 
the regression analysis. The other three heights were equally spaced 
between the maximum and minimum. Rarely was there a significant 
difference (plus or minus 4 meters per second) between the predictions 
calculated by the three orders of polynomial regression (Table XIV). 
Mean velocities predicted using the the third-order regression analysis 
were in error for stations 2S and 5NB at the 13-meter depth. The 
cause of this error could not be explained. 

(3) Confidence Intervals. Confidence intervals (Table 
XV) were computed for the predicted velocities using t values for 
90-percent confidence. In general terms, the confidence intervals were 
larger for all stations in the 1- and 3-meter depths, with the intervals 


becoming less with increasing depth. This corresponds with the lower 
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TABLE XIV 


Predicted Velocities--Depth 
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correlation coefficients found in the shallow depths. Larger confidence 
intervals were also found for data acquired at stations 2N, 2S, 4SA, 
and 4SB. These stations had lower correlation coefficients and higher 
residuals when compared to the other stations. For stations 2N and 2S, 
this may indicate less stability due to the close proximity of the open 
ocean. As stated before, data for stations 4SA and 4SB were acquired 
during spring tide conditions, which tends to indicate less stability for 
prediction. 

In most cases, increasing the degree of the 
regression analysis did not produce more accurate results in terms of 
confidence intervals. The standard error of estimate was sometimes 
lower with a higher order analysis, but the reduction in the degrees of 
freedom countered the effect of the lower standard error of estimate for 
computation of confidence Ince 

Confidence interval results indicate that correcting 
soundings for mean velocity variation based on predictions provided by 
the regression analysis cannot be considered a complete success. 
Attaining ‘the goal of predictions to within plus or minus 4 meters per 
second may not be possible in some cases, but the regression analysis 
does show a significant reduction in error when compared to only using 


average velocity and horizontal variation as the correction procedure. 


B. NOS BAR CHECK DATA 


Analysis of the bar check data revealed too many inconsistencies 
and required too much speculation for: comparison of the CREDDP data 
results. A major problem of the direct comparison method of the bar 
check is that examination of the sound velocity variables and any 
change in their magnitude is not possible. The bar check data did 
indicate the possibility of either an extremely high sound velocity 
variation or an instrument error which was not constant with depth. 
Most likely, the results show an undeterminable combination of sound 
velocity variation and instrument error. The problem of using the bar 


check data for any meaningful analysis is illustrated by Table XVI. 
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TABLE XVI 
Bar Check Correction Data for Survey H-8422 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 


Investigation of sound velocity variation in the Columbia River 


estuary revealed the following significant findings: 


Sound velocity variation in certain areas required the examination 
of how veloctiy changed with _time and location if NOS accuracy 
standards were to be met. The procedure of computing only one 
corrector table based on an average velocity, or correction, was 
not acceptable. 


CREDDP data analysis showed a strong relationship between tide 
eight and sound velocity, which stimulated the development of a 
method to account for the variation of sound velocity with time. 
Using a regression SNE UST the accuracy of VE OCU BnColmeacos 
was Improved significantly, although there were still instances 
where NOS standards could not be met. 


NOS bar check data analysis showed a significant variation. in 
corrections al are the variation did not correlate with tide 
height as the CREDDP data did. Analysis of the variation found 
at each 5S-foot increment of depth revealed inconsistencies which 
could not be explained. 


B. RECOMMENDATIONS 


The following recommendations are made: 


Further See y of sound velocity variation must be carried out for 
verification of ideas and methods presented in this thesis. 
Analysis of data acquired over longer time periods at the same 
locations will allow _a more complete study of the variation . 
associated with different lunar phases and the effect of changes in 
river flow. NOS CTD data were acquired in_the Columbia River 
estuary during the same time period as CREDDP; these data may 
eS UCY of the sound velocity variation over longer time 
periods. 


NOS bar check data analysis showed that the hydrographer should 
be equipped with an efficient method of measuring sound velocity 
in addition to determining instrument error. SDMA the value of 
sound velocity, detection of significant changes and trends is more 
efficient and accurate. Bar checks should not be used as the 
rimary method of correcting soundings, but should be reserved 
or determination of instrument error. 


3) 


Hydrographers conducting surveys in areas where reversing tides 
could cause a significant variation in sound velocity may benefit by 
noting the height of tide for each cast or bar check acquired for 
determining sound velocity corrections. If significant. variation is 
detected, the possibilty of predicting and correcting for temporal 
variation using the tide Hein as an independent variable should 
be investigated. FORTRAN computer programs written for the 
incense: presented in this thesis are listed in Appendices A, B, C, 
, , an 


In areas where temporal variation of sound velocity is. significant, 
the hydrographer ve benefit by using a sound velocimeter to 
correct soundings. he velocimeter would require only a few. 
minutes to deploy and retrieve. The nyene Shap mey could determine 
sound. Vee periodically throughout the working day without 
significantly hindering operations. 
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